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ABSTRACT 

We consider the possibility of discriminating different theories of gravity us- 
ing a recently proposed gravitational wave detector of spherical shape. We 
argue that the spin content of different theories can be extracted relating 
the measurements of the excited spheroidal vibrational eigenmodes to the 
Newman-Penrose parameters. The sphere toroidal modes cannot be excited 
by any metric GW and can be thus used as a veto. 



The efforts aimed at the detection of gravitational waves (GW) started 
more than a quarter of century ago and have been, up to now, unsuccessful 
[H, Resonant bars have proved their reliability, being capable of continous 
data gathering for long periods of time 0, ^. Their energy sensitivity has 
improved of more than four orders of magnitude since Weber's pioneering 
experiment. But a further improvement is still necessary to achieve success- 
ful detection. While further developments of bar detectors are under way, 
two new generations of earth based experiments have been proposed. While 
detectors based on large laser interferometers are already under construction 
P, resonant detectors of spherical shape are under study 0. In the present 
letter we report on a study about the physical information that can be ob- 
tained thanks to the spherical symmetry of the latter detectors. More in 
detail, we want to show that the measurements of the sphere vibrations can 
provide specific information on the field content of the gravitational theory 
predicting the observed features of the waves. 

We would like to remind the reader of the very special position of Ein- 
stein's general relativity (GR) among the possible gravitational theories. 
Theories of gravitation, in fact, can be divided into two families: metric 
and non- metric theories [0. The former can be defined to be all theories 
obeying the following three postulates: 

• spacetime is endowed with a metric; 

• the world lines of test particles are geodesic of the above mentioned 
metric; 

• in local free-falling frames, the non-gravitational laws of physics are 
those of special relativity. 

It is an obvious consequence of these postulates that a metric theory obeys 
the principle of equivalence. More succintly a theory is said to be metric if 
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the action of gravitation on the matter sector is due exclusively to the metric 
tensor. GR is the most famous example of a metric theory. Kaluza-Klein 
type theories, also belong to this class along with the Brans-Dicke theory. 
Different representatives of this class differ for their equations of motion 
which in turn can be deduced from a lagrangian principle. Up to now non- 
Einsteinian and non-metric theories ^ have been considered a curiosity by 
the majority of physicists since there seems to be no compelling experimental 
or theoretical reason to justify their introduction. This point should perhaps 
be reconsidered if we think of the unique role of string theories in this context. 

String theory seems to lead to a consistent framework in which to quantize 
gravity and the other fundamental interactions. Gravity emerges from string 
theory in different forms depending on the chosen vacuum configuration. 
Allowed string vacua are solutions of the equations of motion resulting from 
conformal invariance of the sigma model which governs the string propagation 
[^]. The string equations of motion can also be obtained a posteriori via a 
variational principle. Indeed, the elementary excitations of the string with 
zero mass (light sector) can be described by an effective lagrangian which 
is obtained after integrating out the massive (heavy sector) modes of the 
theory in a procedure a la Wilson. A great variety of string solutions has 
appeared in literature: plane waves, solitons, instantons, rotating and non- 
rotating black holes. More recently a cosmological solution has appeared that 
seems to lead to a potentially detectable quantity of gravitational radiation 
in a wide spectrum of frequencies 0. These solutions can be both metric 
and non-metric but are certainly non-Einsteinian since all of these vacua 
contain a massless scalar called dilaton. When perturbations of these vacua 
are considered all fields fluctuate and the resulting theory is non-metric. But, 
as in string theory the expectation value of the dilaton plays the role of a 
coupling constant, it is generally believed that some non-perturbative effect 
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(the most serious candidate at the moment is supersymmetry breaking) gives 
a mass to the dilaton. If this is true, the dilaton interactions with the other 
fields may well be negligible [0 and the theory may revert to be metric. 



These are our motivations for testing theories of gravity using the mea- 
surable toroidal and spheroidal vibrational eigenmodes of a 

sphere. The signature of a non-metric theory could obviously also be 
detected in experiments on the equivalence principle. 

Before discussing the interaction with an external GW field, let us con- 
sider the basic equations governing the free vibrations of a perfectly homo- 
geneous, isotropic sphere of radius i?, made of a material having density p 
and Lame coefficients A and /i ITT . 



Following the notation of [jl2||, let Xj, i = 1, 2, 3 be the equilibrium position 
of the element of the elastic sphere and x'^ be the deformed position then 
the displacement vector. Such vector is assumed small, 
so that the linear theory of elasticity is applicable. The strain tensor is 
defined as Uij = (1/2) (ujj- + Uj^i) and is related to the stress tensor by aij = 
6ijXuii + 2puij. The equations of motion of the free vibrating sphere are thus 



with the boundary condition: 



{SijXuii + 2/iMij) (1) 



nj(Ty = (2) 

at r = i? where = Xi/r is the unit normal. These conditions simply 
state that the surface of the sphere is free to vibrate. The displacement 
Ui is a time-dependent vector, whose time dependence can be factorised as 
Ui{x,t) = Ui{x)exp{iujt), where uj is the frequency. The equations of motion 
then become: 

/iV^Ui + (A + p)'Vi{'VjUj) = -uP'pui (3) 
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Their solutions can be expressed as a sum of a longitudinal and two transverse 
vectors: 

u{x) = CoV0(f) + CiLxix) + C2V X Lx{x) (4) 

where Co,Ci,C2 are constants and L = x x V is the angular momentum 
operator. Regularity at r = restricts the scalar functions and x to 
be expressed as (l){r,9,ip) = ji{qr)Yim{0, ip) and x{r,0,ip) = ji{kr)Yi^{e,ip). 
Yim{0,ip) are the spherical harmonics and ji the spherical Bessel functions 

[0: 

. , , /I d X'- /sinx'" 
\xdxj V X 

= puj"^ / (A + 2yu) and k"^ = puj"^ / p are the longitudinal and transverse wave 
vectors respectively. 

Imposing the boundary conditions (0) at r = i? yields two families of 
solutions: 

• Toroidal modes: these are obtained by setting Cq = C2 = 0, and 
Ci 7^ 0. The eigenfunctions have the form: 

ulUr.O,^) = T^i{r)LYUe,^) (6) 

with Tni{r) proportional to ji{knir). The eigenfrequencies are deter- 
mined by the boundary conditions (0) which read 

h{kR) = (7) 

where 



dz 



(8) 



z 

Spheroidal modes: these are obtained by setting Ci = 0, Cq 7^ and 
C2 7^ 0. The eigenfunctions can be conveniently rewritten as 

Unirni.^) = Aniir)Yimi0, v)^ " Bni{r)n X LYi^{e, (9) 
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where Ani{r) and Bni{r) are dimensionless radial eigenfunctions [|T2 
which can be expressed in terms of the spherical Bessel functions and 
their derivatives. The eigenfrequencies are determined by the boundary 
conditions which read 

fUqR) - l^q^R'UqR) l{l + l)fi{kR) \ ^ ^ 

' V fi{qR) '^MkR) + - l]fo{kR) ) 

(10) 

where 



/0(^) = ^ (11) 



The eigenfrequencies can be determined numerically for both toroidal and 
spheroidal vibrations. Each mode of order / is (2/ + l)-fold degenerate. In 
the table below we show the value of the {kR) roots for the lowest toroidal 
and spheroidal modes of vibration. 



/ 


n 


{kR) toroidal 


( kR) spheroidal 


I 


n 


{kR) toroidal 


{kR) spheroidal 





1 




5.4322 


3 


1 


3.8647 


3.9489 




2 




12.138 




2 


8.4449 


6.6959 




3 




18.492 




3 


11.882 


9.9720 




4 




24.785 




4 


15.175 


12.900 


1 


1 


5.7635 


3.5895 


4 


1 


5.0946 


5.0662 




2 


9.0950 


7.2306 




2 


9.7125 


8.2994 




3 


12.323 


8.4906 




3 


13.211 


11.324 




4 


15.515 


10.728 




4 


16.544 


14.467 


2 


1 


2.5011 


2.6497 


5 


1 


6.2658 


6.1118 




2 


7.1360 


5.0878 




2 


10.951 


9.8529 




3 


10.515 


8.6168 




3 


14.511 


12.686 




4 


13.772 


10.917 




4 


17.886 


15.879 
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The eigenfrequency values can be obtained from : 



i^nl — 




1^ {kR)nl 

p R 



(12) 



The detector is assumed to be non-relativistic (with sound velocity f s -C c 
and radius i? -C A the GW wavelength) and endowed with a high quality 

factor {Qni = i^niTni ^ Ij whcrc Tni is the decay time of the mode nl). The 
displacement m of a point in the detector can be decomposed in normal modes 
as: 



N 

where N collectively denotes the set of quantum numbers identifying the 
mode. The basic equation governing the response of the detector is 



We assume that the gravitational interaction obeys the principle of equiva- 
lence which has been experimentally supported to high accuracy. In terms 
of the so-called electric components of the Riemann tensor E^j = Roioj, the 
driving force fN{t) is then given by 



where M is the sphere mass and we consider the density p as a constant. In 
any metric theory of gravity Eij is a 3 x 3 symmetric tensor, which depends 
on time, but not on spatial components. 

Let us now investigate which sphere eigenmodes can be excited by a 
metric GW, i.e. which sets of quantum numbers N give a non-zero driving 
force. 

a) Toroidal modes 




(13) 



(14) 




(15) 
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The displacement vector can be expressed as in eq. (|^). Up to a normalisation 
constant C, the driving force is 



E,n, — E T. 



drr^ji[k^\) / d6 sin 6 



^yy 



sin ^ sin 20^^ + cos ^ cos 20 '™ 



89 

BY* dY* 
+ E^y [ sin 9 cos 20^:^ - cos 9 sin 20 



90 



+ Exz 

+ Eyz 

+ (e. 



— sin cos 9 



89 
8Y, 



Im 



89 



+ (sin 6* cos ( 



cos^ 9 
sin 9 



8Y, 



COS( 



Im 



COS COS 6'- 



+ (sm sm — - 

ao' sm 



Im 



8Y 



sm I 



yy 



cos9^ 



Using the equations 



89 



2l + l{l-m)\ 



An {I + m) 



l 8Pr{cos9) ^_^^, 



89 



and 



8Y, 



Im 



-im[ 



21 + 1 (/ -m)!]5 



P;'"(cos^)e-^'"^ 



(16) 



(17) 



(18) 



An {l + m)\. 

the integration over can be performed. Eq. ([T6| ) then contains integrals 
over 9 of the form: 

.8P^Hcos( 



[8111^9 -cos"^ 9) P[^\cos( 



sin 6' cos 6*- 



89 



d9 



and 



2^9P±2(C0S 



2 sin cos 9 PPi cos 9) + sin^ 9 



d9 



(19) 



(20) 



' ' 89 

After integration by parts, the derivative terms in eqs. ([T9|) and (^) exactly 
cancel the non-derivative ones. The remaining boundary terms vanish too, 
thanks to the periodicity of the trigonometric functions and to the regularity 
of the associated Legendre polynomials. The vanishing of the above integrals 
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has a profound physical consequence. It means that in any metric theory of 
gravity the toroidal modes of the sphere cannot be excited by GW and can 
thus be used as a veto in the detection. In this respect we disagree with the 
results of that seems to find a zero result only for the even / case. 



b) Spheroidal modes 
The forcing term is given by: 

(21) 

One is thus lead to compute integrals of the following types 

J x^x%miO,^)d^x (22) 

and 

j x^x'LuYi^{e,ip)(fx (23) 

Since the product x'^x^ can be expressed in terms of the spherical harmonics 
with Z = 0, 2 and the angular momentum operator does not change the value 
of /, one immediately concludes that in any metric theory of gravity only the 
/ = 0, 2 spheroidal modes of the sphere can be excited. At the lowest level 
there are a total of five plus one independent spheroidal modes that can be 
used for GW detection and study. 



These results agree with the conclusion of |0 and generalize the results 



of PI, |T6 



From the analysis of the spheroidal modes active for metric GW, we now 
want to infer the field content of the theory. For this purpose it is convenient 
to express the Riemann tensor in a null (Newman- Penrose) tetrad basis P]. 

To lowest non-trivial order in the perturbation the six indipendent "elec- 
tric" components of the Riemann tensor may be expressed in terms of the 
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Newmann-Penrose (NP) parameters as 



Im^4 Re^4 - $2: 



22 



Jm\l/4 



22 



-2v/2i?e^3 \ 

2^2/771^3 

-6^2 / 



(24) 



The NP parameters allow the identification of the spin content of the metric 
theory responsible for the generation of the wave 0. The classification can 
be summarized in order of incresing complexity as follows: 

• General Relativity (spin 2): \E'4 7^ while ^^2 = ^E's = $22 = 0. 

• Tensor-scalar theories (spin 2 and 0): \I'4 7^ 0, \E'3 = 0, \I'2 7^ and/or 
$22 ^ {e.g. Brans-Dicke theory with \&4 7^ 0, \&2 = 0, \I^3 = and 
$22 7^ 0). 

• Tensor- vector theories (spin 2 and 1): \I'4 7^ 0, \E'3 7^ 0, $22 = ^'2 = 0. 

• Most General Metric Theory (spin 2, 1 and 0): ^4 7^ 0, ^^2 7^ 0, ^3 7^ 
and $22 7^ 0, {e.g. Kaluza-Klein theories with ^^4 7^ 0, \E'3 7^ 0, $22 7^ 
while \E'2 = 0). 

In eq. (0), we have assumed that the wave comes from a localized 
source with wave vector k parallel to the z axis of the detector frame. In 
this case the NP parameters (and thus the wave polarisation states) can be 
uniquely determined by monitoring the six lowest spheroidal modes. If the 
direction of the incoming wave is not known two more unknowns appear in 
the problem, i. e. the two angles of rotation of the detector frame needed to 
align k along the z direction. In order to dispose of this problem one can 
envisage the possibility of combining the pieces of information from an array 
of detectors [|17|. We restrict our attention to the simplest case in which the 
source direction is known. 
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The NP parameters of the incident wave can be easily obtained from the 
following argument. It is well known that any 3x3 symmetric tensor such 
as (|2^) can be decomposed in the following way 



l,m 



{l,m) 



(25) 



(2,m) 



where Slj' = 5ij/y/ATT (with 5ij the Kronecker symbol) and 5*, 
—2, ..2) are five linearly independent symmetric and traceless matrices. The 
following explicit representation is particularly suited for exposing the spin 



content of the GW \12 



S. 



(0,0) 
^3 



S. 



(2,2) 



s. 



(2,1) 



, 15 , 1 

Y 

^327r^ 

-(^ 

^327r 



/I 0\ 

1 
lo 1/ 

/I i 0\ 
i -1 

Vo 0/ 

/O 1\ 

)^ z 

VI ^ 0/ 



(2,0) 



IGtt' 



s. 



(2,-2) 



5, 



(2,-1) 



, 15 , 
^327r' 



, 15 . 
^327r' 



/-I 
-1 

V 
/ 1 











2/ 

— i 

i -1 
V 0/ 

/O 1 \ 

-z 
\l -I / 



(26) 



The matrices S'^ are trace-orthogonal and are connected to the spherical 
harmonics by 





— ^0,0? 


q(2,0) 


- ^2,0 


r-(2,2) 


= ^2,2, 


r-(2,-2) 


= ^2,-2 


r.(2,l) 


= >^2,1, 


r.(2,-l) 


= 



(27) 



The vector in eqs. (^) has been defined after eq. 
Taking the scalar product we find 

T 



Co,o(i) 



-SifEij{t) 
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C2,m{t) = —Sij'^Eij{t) 

For the NP parameters we find 

Im^i = - 

V 

JmWs = — [C2,l + C2 -1 



15 



(28) 



15 



C2.2 + C2 -2 



(29) 



V 327r^ ' ' " V 327r' 

1/15. 1 . T ^ 

leVvt^^'^-^^'-^] '"^^^=16\ n 

Eqs. (|2^) relate the quantities Q^m with the GW polarization states, de- 
scribed by the NP parameters. A pictorial representation of the six po- 
larization states connected with the NP parameters is given in Figure 1 
where the wave vector is assumed parallel to the z axis as in (|2^). Start- 
ing from the top left corner of the figure we draw the polarization states 
^2, $22, -Re^4, /m^4, -Re^^s, 




Figure 1 
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Eq. (|29| ) can be put in correspondence with the output of experimental 
measurements if the ci^m are substituted with their Fourier components at 
the quadrupole and monopole resonant frequencies which, for the sake of 
simphcity, we collectively denote by ujq. The Q^m(i^o) can be determined in 
the following way: once the Fourier amplitudes An{ujo) are measured, by 
Fourier transforming (|T^) and (|T^) we get the Riemann amplitudes Eij{uJo) 



which, using (|28D , yield the desired result. 

In order to determine the An^ujq) amplitudes from a given GW signal 
two conditions must be fulfilled: 

• the vibrational states of the five-fold degenerate quadrupole and monopole 
modes must be determined. The quadrupole modes can be studied by 
properly combining the outputs of a set of at least five motion sen- 
sors placed in independent positions on the sphere surface. Explicit 
formulas for practical and elegant configurations of the motion sensors 
have been reported by various authors [0, |TB], |Ty|. The vibrational 
state of the monopole mode is provided directly by the output of any 
of the above mentioned motion sensors. If resonant motion sensors are 
used, since the quadrupole and monopole states resonate at different 
frequencies, a sixth sensor is needed. 

• The spectrum of the GW signal must be sufficiently broadband to over- 
lap with the antenna quadrupole and monopole frequencies. 

From table 1 we see that the second order quadrupole spheroidal mode is 
close to the lowest order monopole and also to the toroidal mode n = 1, 



/ = 4. Since it has been demonstrated ||20| that the second order quadrupole 



spheroidal mode has high cross section for general relativistic GW (only a 
factor 2.6 lower than the cross section of the first order quadrupole spheroidal 
mode), the full analysis of the GW signal for the most general metric GW, 
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including the toroidal mode veto, can be hopefully performed in a small fre- 
quency range. The above reported analysis can then be applied to a large 
class of GW sources, including gravitational collapses, inspiralling and coa- 
lescing binary systems and stochastic sources. 

The sensitivity of a multiton spherical detector, making use of ultracryo- 
genic and superconducting techniques for noise reduction could be such to 



detect these events with a reasonable rate ||14 |. 
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